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Radio flares and plasmon size in Cygnus X-3 
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ABSTRACT 

We have observed a number of minor radio flares in Cyg X-3 using the MERLIN array. 
Photometric observations show the system to be highly active with multiple flares on 
hourly timescales over the one month observing programme. Analysis of the source's 
power spectrum show no persistent periodicities in these data, and no evidence of the 
4.8 hr orbital period. An upper limit of 15 mJy can be placed on the amplitude of any 
sinusoidal varaition of source flux at the orbital period. The brightness temperature 
of a flare is typically Tb > 10 9 -10 10 K, with a number of small flares of 5 minute 
duration having brightness temperatures of Tb > few x 10 11 K. For such a change in 
flux to occur within a typical 10 minute timescale, the radiation must originate from 
plasmons with a size < 1.22 AU. This emission is unlikely to originate close to the 
centre of the system as both the jets and compact object are buried deep within an 
optically thick stellar wind. Assuming a spherically symmetric wind, plasmons would 
become visible at distances ~ 13 AU from the core. 

Key words: binaries: close, stars: individual: Cyg X-3, stars: Wolf-Rayet, stars: 
variables: other, radio continuum: stars, X-rays: stars 



1 INTRODUCTION 

The energetic X-ray binary Cygnus X-3 has been studied in 
great detail since its discovery in 1967 (Giacconi et al. 1967). 
It is an active source with measurements in the radio, sub- 
mm, infrared, X-ray and gamma-rays all of which indicate 
a highly variable source. 

Following infrared spectroscopic measurements (van 
Kerkwijk et al. 1992; Fender, Hanson & Pooley 1999), the 
accepted morphology of Cyg X-3 is a compact object in or- 
bit around a Wolf-Rayet star of the WN4-5 type. A mea- 
sured infrared and X-ray orbital period for the binary of 
4.8 hr indicates a very compact system (< 10 Rq) within 
which a large degree of interaction takes place. The proxim- 
ity of a compact object to a star with a strong wind provides 
both accretion onto the compact object, but also a medium 
of varying optical depth through which jets, formed from 
accreted material, must travel. Any emission from the jets 
cannot be observed directly until the wind becomes optically 
thin to the radiation in the jets. This creates a time-lag be- 
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tween radiation at different frequencies in the jet (typically 
radio and sub-mm emission) and between jet and central 
infrared and X-ray emission. 

One can use the different radiation regimes to inves- 
tigate different properties of the system. For example, the 
low-frequency radiation at radio and sub-mm frequencies 
comes from synchrotron emitting electrons in the jets, and 
this has been used by Ogley et al. (in prep.) to obtain an 
upper-limit to the magnetic field in the jets of 260 /iT (2.60 
G) at 100 Rq from the core. As one observes higher fre- 
quencies, the emission changes from synchrotron to thermal 
free- free, which now comes from the wind in the system and 
not the jets. An example of the wind emission observed by 
the ISO satellite is given in Ogley, Bell Burnell & Fender 
(2001), and a complete spectrum showing the different emis- 
sion mechanisms from radio to infrared is given in Ogley et 
al. (in prep.) Fig. 6. However, a thorough investigation us- 
ing a single radiation regime can also provide a great deal 
of information. 

The accretion process is not constant, but is probably 
variable in both density and mass-transfer rate. This causes 
a change in the number of electrons which are available for 
acceleration in the jet, which alters the synchrotron luminos- 
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ity. Radio photometry monitors the change in jet emission 
following a change in the Wolf-Rayet state and three dis- 
tinct states of emission in the unresolved core of the system 
have been identified (Waltman et al. 1995): 

(i) For the majority of the time the system is in a quies- 
cent radio state with GHz flux densities around 100 mjy. 

(ii) Once or twice a year the system undergoes a major 
flare. The radio flux density increases by 10-100 times on a 
timescale of around a day. A major flare is very energetic 
and lasts for days or weeks, exponentially decaying back to 
quiescence. Identified immediately before a major flare is a 
period of quenching during which the flux density drops to 
around 10 mjy and remains at this lower value for up to a 
week. 

(iii) An intermediate state exists between the high flux 
of a major flare ( > 1 Jy) and the normal quiescent flux of 
100 mjy. This is the state of minor flaring and can last for 
several weeks. During a minor flare state, the system can 
undergo multiple short-duration flares and is highly active. 

Continuous monitoring determines the radio state of 
Cyg X-3, and allows detailed investigation of the jet mech- 
anism and structure during these different states. 
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Figure 1. 15 GHz radio photometry from the Ryle telescope 
which triggered the start of our MERLIN observation programme. 
A minor flare occurred on MJD 50407 and our MERLIN obser- 
vations were triggered from MJD 50419 until MJD 50459 as dis- 
cussed in the text. See also Fig. 2. 



Table 1. Date and calibration sources for the MERLIN observa- 
tions. 



2 THE RADIO OBSERVATIONS 

Daily radio monitoring was carried out using the Ryle Tele- 
scope in Cambridge and the Green Bank Interferometer 
(GBI) in West Virginia. These observatories provided pho- 
tometry at 2.3 and 8.3 GHz (GBI) and 15 GHz (Ryle). Iden- 
tification of a minor flare in 1996 November triggered a set of 
interferometric observations using the MERLIN array. The 
MERLIN observations were scheduled at C-band (5 GHz) 
with 6 epochs of 12 hours in duration each separated by one 
week. The aim of these observations was to map the emis- 
sion from a plasmon as it travelled from the centre of the 
binary, as observed by Newell (1995). At a frequency of 5 
GHz a plasmon at 10 kpc, travelling at a speed of 0.35 c in 
the plane of the sky, would separate from the core at one 
MERLIN beam- width per epoch (~ 50 mas). This jet veloc- 
ity, taken from the apparent transverse motion measured by 
Spencer et al. (1986) represents a lower-limit to the separa- 
tion per epoch. Other measurements have indicated higher 
velocities of components close to the plane of the sky (~ 0.6c 
from Marti' et al. (2000) and ~ 0.8c from Mioduszewski et 
al. (2000)). These larger apparent velocities would therefore 
yield greater separations per MERLIN epoch. 

Fig. [j] shows the radio photometry from the Ryle tele- 
scope at 15 GHz which was used as a trigger for the MERLIN 
observations. Cyg X-3 was in quiescence until MJDp] 50400 
when it underwent a mild quenching period (flux around 50 
mjy). Shortly after the quenching, a minor flare occurred on 
MJD 50407 which signalled the start of a minor flare period. 
Subsequent minor flares occurred at random times after this 
date. 

Observations were taken using the standard 5 GHz 

T In this paper we use the definition of modified Julian date to be 
MJD = JD - 2 400 000.5, so integer values of MJD fall on 0000 
UT. 



Calibrators 

MJD Date Point source, flux Phase 



50419 


01 December 1996 


OQ 208 


2005+403 


50425 


07 December 1996 


OQ 208 


2005+403 


50432 


15 December 1996 


0552+398 


2005+403 


50438 


21 December 1996 


OQ 208 


2005+403 


50446 


28 December 1996 


OQ 208 


2005+403 


50459 


11 January 1997 


0552+398 


2005+403 



MERLIN continuum setup. A 15 MHz bandwidth was used 
for observations of the target source Cygnus X-3, a phase 
calibrator source 2005+403, a flux calibrator 3C84, and a 
point source calibrator, either OQ 208 or 0552+398. To cal- 
ibrate the flux easily, observation of a point source of known 
flux is required. However, on the longer MERLIN baselines 
all the point sources are variable in flux, and all the constant 
flux sources are resolved. One has to observe a point source 
over all baselines and calibrate the flux scale by measure- 
ments of a source of known flux using the short baselines 
only. Details of this calibration procedure are given in the 
MERLIN user's guide, Thomasson et al. (1993) and Ogley 
(1998). 

The source and observational details for all the epochs 
are given in Table hi The times of the MERLIN observations 
are shown in Fig. g together with Ryle photometry at 15 
GHz around the same time. The MERLIN observations are 
indicated by the arrows. 



3 MERLIN PHOTOMETRY 

Because of the variable nature of the source at the time 
of the MERLIN observations, a robust map of the source 
is almost impossible to obtain. Variations in observed am- 
plitude can be due to two things: a variable core; or emis- 
sion away from the phase centre. To reconstruct an image, 
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Figure 2. Times of the MERLIN observations indicated by ar- 
rows overlaying a plot of Ryle photometry at 15 GHz. 



the assumption that any variation in the fringe amplitude 
with hour angle is due to the interferometer response to the 
source structure only is made. However, this is violated if 
the source varies within the time needed for aperture syn- 
thesis. For our observations, a direct reconstruction of the 
source results in many artifacts including jet-like emission 
from a central core. The variability of the source is shown in 
Fig. ^ which plots the measured flux against baseline length 
or uv spacing. One can see that the source is highly variable 
over all baseline distances and no accurate reconstruction of 
the source brightness distribution can be obtained. 

The plot shown in Fig. ^ does allow some quantative 
evaluation. A resolved source would show a decrease in mea- 
sured flux with increasing uv spacing, whereas during a flare 
the measured flux is approximately constant over all uv sep- 
arations. In Fig. 3 the measured flux is high at the longest 
baselines suggesting that the source is unresolved. A source 
that is unresolved at a uv spacing of 3.5 MA at 5 GHz has 
an angular size of less than 70 mas. This is consistent with 
a scatter-broadened size of Cygnus X-3, which has a strong 
dependence on frequency (Wilkinson, Narayan & Spencer 
1994) and has been measured at 5 GHz to be ~ 20 — 30 mas 
(Mioduszewski et al. 2000). 

For these observations the success of the MERLIN array 
is not in its imaging capabilities, but its photometry. Figs 
^]|] show photometry for all six C-band epochs. Note that 
in these figures the vertical scales vary and the zero is not 
included. Vertical lines in all these figures indicate the time 
at which X-ray minimum occurs using the ephermeris by 
Matz et al. (1996). Data have been averaged into time bins 
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Figure 3. A plot of flux against uv distance for the entire first 
epoch. There is significant, and unstructured, variation of flux 
which is unlikely to be caused by source structure and is due to 
changes in the source flux over the epoch. 



Hours from MJD 50419.0 
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Figure 4. Cyg X-3 photometry on 01 December 1996. Time bins 
are 3.4 min. Vertical lines indicate time of X-ray minimum from 
the ephermeris by Matz et al. (1996). 
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Figure 5. Cyg X-3 photometry on 07 December 1996. Time bins 
are 2.2 min. 
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Figure 6. Cyg X-3 photometry on 15 December 1996. Time bins 
are 3.8 min. 
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Figure 7. Cyg X-3 photometry on 21 December 1996. Time bins 
are 4.0 min. 
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Figure 8. Cyg X-3 photometry on 28 December 1996. Time bins 
are 2.0 min. 
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Figure 9. Cyg X-3 photometry on 11 January 1997. Time bins 
are 3.4 min. 
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frequency [l/days] 

Figure 10. A raw power spectrum from all six epoch photomet- 
ric points. The insert is the Fourier transform of the sampling 
function of the data. 

of 2.0-4.0 minute duration to reduce noise and all baselines 
have been used. 

A number of features are visible in the plots. Points to 
note in the data are: 

(i) In all but the first four hours of the MJD 50445 (28 De- 
cember 1996) observation, the flux is variable on timescales 
shorter than a couple of hours. This causes problems in the 
mapping as discussed above. 

(ii) A number of flares occurred, however, the timescales 
involved are not constant or correlated with flare strength 
or length of time in quiescence preceding a flare. 

3.1 Power spectra 

To search for any periodic signals in the data, photometric 
points from all six epochs were included and the Fourier 
spectrum was determined. This is shown in Fig. [l(]. The 
main part of the figure is a direct Fourier transform of the 
data. The insert is the Fourier transform of the windowing 
function. 

The main transform of the data is the time-domain 



© 0000 RAS, MNRAS 000, 000-000 



Radio flares and plasmon size in Cygnus X-3 5 



4.1 Brightness temperature of the flares 




frequency [l/'days] 



Figure 11. A cleaned spectrum. The insert shows the power at 
a period of 4.8 hours. There is no evidence for any periodic signal 
at this binary period. 




frequency [l/days] 



Figure 12. Power spectra of the data cleaned in an identical way 
to Fig. but split into two halves. The top panel shows the first 
three epochs, and the bottom panel shows the latter three epochs. 
No correlation in the periods are found between the two halves 
of the data. 



equivalent of an interferometric dirty map, and the trans- 
form of the windowing function is the equivalent of a dirty 
beam. In the same way as one can re-create a clean map 
from a dirty map and dirty beam, the power spectrum can 
be cleaned by an iterative subtraction of the window func- 
tion (Clark 1980; Roberts, Lehar & Dreher 1987). The result 
of this procedure is shown in Fig. [ll]. 

To check whether the periods shown in Fig. [ll] are spu- 
rious, the first three epochs were analysed separately from 
the last three epochs. The clean spectra are shown in Fig. 
|li| . As there is no correlation between the two halves of the 
data sets, we can conclude that there are no stable period- 
icities and any periodicity shown is due to spurious features 
or quasi-periodic oscillations. 



4 DISCUSSION 



A brightness temperature of a flare can be calculated asso- 
ciated with the flux change from the source. A black body 
sphere with radius r, flux S v , at a distance D, assuming a 
Rayleigh- Jeans tail, would have a brightness temperature of 



2 S V D 2 



(1) 



where c is the speed of light, fee is Boltzmann's constant and 
v is the frequency. Therefore a flare of duration At minutes 
and flux change of AS mjy has a brightness temperature of 



T b = 9.66 x 10 1 



ASmJyD, 



mjy -f^kpc 

A/;' : . 



K 



(2) 



where quantities are measured in their subscripted units. 

For a typical flare: MJD 50419.56, the flux increase 
AS = 80 mjy over At = 46 minutes at v = 5 GHz. If Cyg 
X-3 is at a distance of 10 kpc, the brightness temperature 
is Tb > 1.5 x 10 10 K. This value is typical for the major- 
ity of flares, including the large flare at MJD 50432. The 
largest brightness temperature occurred at MJD 50438.6, if 
the change within 1 integration bin is believable. For this 
flare the flux change is AS = 10 mjy over At = 4.0 min, 
and so the brightness temperature is Tb > 2 x 10 11 K. 

There is a maximum brightness temperature above 
which inverse Compton losses become catastrophic. This 
maximum brightness temperature can be written for galac- 
tic sources as 



< 1.6 x 10 1 



^ 0.8 


Hi -a) 


vv~ 








V 





V 1 



(3) 



where vu is the upper frequency of the synchrotron emission 
in GHz, T> is the Doppler boosting factor of the emission 
and c^g' is the a-dependent parts of the ratio of Pacholczyk 
constants C5 and ce where 

3a + 5 r 



„(«) 



3a + 1 

6 



6q 2 + 19a + 13 r ( 6a + 5 ) r ( 6a + 13 ) 



(4) 



(Pacholczyk 1970; Hughes & Miller 1991). For a spectral in- 
dex of a = —0.1, and vu = 350 GHz (Ogley et al. in prep.), 
and a frequency of v = 5 GHz used in this work, the ex- 
pression for the maximum brightness temperature evaluates 
to 

T max < 4.8 x lO 11 © 12 K. (5) 

An estimation of the bulk Doppler factor of the jets can be 
inferred from jet expansion. An arcsec-scale measurement 
(Marti et al. 2000) has shown that at that distance from 
the core, components are expanding with proper motions of 
fi a = 10.9 mas d _1 and fi r = 9.8 mas d _1 for the approach- 
ing and receding components respectively. This allows an 
estimate of the jet bulk velocity and angle to the line of 
sight to be measured, where 



/3cos( 



(6) 



Ma + fir 

From these observations the authors find /3cos# = 0.053 ± 
0.027, from which they assume a jet velocity of j3 = 0.60 and 
an angle to the line of sight of 8 = 85°. If these values are 
adopted then a Doppler factor of 27 = 1.18 can be calculated. 
If the jet speed is increased to /3 = 0.89 as observed in VLBA 
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observations (Mioduszewski et al. 2000), then by using the 
same value of 0, the Doppler factor increase to T> = 2.02. 
Using this range of Doppler factors, the maximum brightness 
temperature is therefore T max < 6 — 11 x 10 11 K. 

4.2 Orbital motion 

The effect of the binary's orbital motion on the observed 
radio emission from Cyg X-3 is found to be negligible, unlike 
Cyg X-l (Brocksopp et al 1999). Reasons why this could be 
are 

(i) The binary is face on (which is irrespective of the jet 
angle). Throughout the orbit there is no difference in the 
amount of material that the jets encounter which would scat- 
ter or absorb radiation. 

(ii) Radio emission in Cyg X-3 is dominated by optically 
thin events which occur from a region way outside the or- 
bital diameter, and so any 4.8 h variation, in say particle 
production or blob ejection, is smoothed out by the time we 
see it. In Cyg X-l, the radio emission is self-absorbed and is 
dominated by the high frequencies, this occurring closer to 
the binary centre. 

To quantify this absence of radio emission, each run was 
Fourier transformed into 100 frequencies with a separation 
of 0.1 cycles per day. The power spectra were then summed 
over all the data sets. No significant peak at or around 5 
cycles per day (= 4.8 h period) was found in either the 
individual spectra or in the sum. A 2 a upper limit of 15 mjy 
can be set for the amplitude of any sinusoidal variation in 
source flux for the average power spectrum. Each daily run 
was also split into 5 overlapping data sets, transformed and 
averaged to give an average power sepctrum per run. These 
showed significantly higher power at low frequencies with 
a power-law dependence with a typical slope of —2.4 ± 0.2 
over the frequency range 4-40 cycles per day. The power 
law shows that no particular timescale can be assigned to 
the data. 

We have indicated on the figures ^-|^, by vertical lines, 
the times of X-ray orbital minimum as calculated from the 
ephemeris of Matz et al. (1996), and one can see that any 
radio state is equally probable. 

4.3 Emission position 

From the small time-scale variations in the photometric data 
flares can be identified as originating from a region typically 
10 light minutes, or 1.22 AU across. However, it is highly 
unlikely that this emission occurs close to the centre of the 
system where the jets are accelerated, as both the jets and 
the compact object are buried deep within an optically thick, 
dense stellar wind from the companion star, see Ogley, Bell 
Burncll, Fender (2001). 

The varying opacity to the radiation from the wind cre- 
ates zones along the jet where radiation at different frequen- 
cies dominates, as the surrounding medium becomes opti- 
cally thin. Radii at which this occurs can be calculated using 
a modified Wright & Barlow (1975) wind model, a function 
for the Gaunt factor appropriate for radio wavelengths (Lei- 
therer & Robert 1991) and assuming spherical symmetry. 
The radius at which emission occurs then depends on only 



the wind parameters such as composition, temperature and 
mass-loss rate, as well as the frequency of the observed emis- 
sion. From the delay of flares at 2.25, 8.3 and 15 GHz during 
1991 to 1995, parameters for the mean number of free elec- 
trons, 7 = 1, the RMS ionic charge, Z — 1, the mean atomic 
weight per nucleon, fi = 4, the temperature, T — 20, 000 and 
the mass-loss rate, M = 1 X 10 -5 Mq yr _1 can be calculated 
(Waltman et al. 1996). This then simplifies the frequency- 
dependent radius to 

R c = 3.28 x 10 3 j 1/3 ^ 3 Rq, (7) 
where 

g = 9.77 + 1.27 log (2.83 x 10 3 v^) (8) 

and z/ghz is the frequency measured in GHz. For a frequency 
of 5 GHz, the radius at which emission is observed is there- 
fore R c > 13 AU. This creates significant layering in the 
emission from the jets with the higher frequencies occurring 
closer to the centre (at 350 GHz, emission would occur at 
only 0.5 AU; see Ogley et al. in prep.). We note that the 
Waltman et al mass loss rate assumes an expansion veloc- 
ity of only 0.3c; a higher velocity would give a larger mass 
loss rate and larger values of R c . Furthermore Ogley et al 
(2001) find a significantly larger mass loss rate from infrared 
observations, a rate independent of the expansion velocity. 
Finally we note that the wind may not be spherically sym- 
metric - see Ogley, Bell Burnell, Fender (2001) and Fender, 
Hanson, Pooley (1999). 



5 CONCLUSIONS 

We have observed the X-ray binary Cygnus X-3 with the 
MERLIN interferometer with the intention of mapping jet 
components in a period of minor flare activity. We observed 
for six epochs between 01 December 1996 and 11 January 
1997. In all epochs the source showed some degree of vari- 
ability and consequently we were unable to map the source 
without the creation of artifacts. 

High time-resolution photometric observations show a 
number of small-flux flares of around an hour in duration, 
with a large amount of structure. This is typical of a minor- 
flare period. Power spectra of the data show no persistent 
periodicities in the data, and no evidence of the orbital pe- 
riod of 4.8 hr. 

A measurement of the brightness temperature for the 
flares show typical values of 10 9 -10 10 K, with the largest 
values occurring over 4 minutes and a brightness tempera- 
ture of 2 x 10 11 K. The flare emission is from a region of 
typically 10 light minutes across, with a diameter of 1.22 
AU. However, to be visible at a frequency of 5 GHz these 
plasmons would be situated at a distance of 13 AU from the 
core, assuming a spherical wind. 
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